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Abstract 



Q_ij The power-law tail of high-pj it spectra observed in forward d+Au collisions at 

RHIC can be attributed to the power-law decrease of the dipole forward scattering 
amplitude appearing in the color glass condensate (CGC) approach. Forward particle 
production probes the small-x gluon distribution of the target nucleus where its 
anomalous dimension is rather flat (7 = 0.6 ~ 0.8) for moderately high pt (^ 5 
GeV), and where the leading-twist DGLAP approximation is not valid. In the same 
framework, we examine p and p production using baryon fragmentation functions 
parameterized in the Lund fragmentation scheme. This provides a good description 
of the forward p spectrum while it underestimates the p data by as much as a factor 
of 2 ~ 3 at pt ;$ 4 GeV. Part of this anomalous baryon excess can be attributed to 
surviving constituent diquarks from the deuteron projectile. Thus, the contribution 
from diquark scattering may play an essential role for forward baryon formation. 



1 Introduction 

One mysterious problem in high-energy p(d)+A collisions is a baryon excess seen in the 
p/ir + ratio, observed at semi-hard transverse momenta (pt^t 1 GeV) over a wide range of 
rapidity. Measurements include the central rapidity region at Fermilab (\fs = 27.4, 38.8 
GeV) |'l] and RHIC (x/s = 200 GeV) as well as forward rapidities (y ~ 3) at RHIC 0. 
In a phenomenological picture, the standard string models assume a diquark-quark string 
structure of the proton beam and dominance of soft processes for leading baryon produc- 
tion as in the inside-outside cascade picture. This leads to baryon production only near 
beam rapidity and with small p t . To explain the baryon production near midrapidity two 
additional nonperturbative models have been proposed: the diquark breaking jlj and gluon 
junction [5j mechanisms. However, these models are also restricted to small momentum 
transfer (p t <, 2 GeV). 

Next-to-leading order (NLO) perturbative QCD calculations with standard fragmenta- 
tion functions cannot be reconciled with the forward-rapidity data from d+Au collisions 
at RHIC, either. The data clearly show a large baryon excess jBJ- Meanwhile, a pointlike 
diquark picture involving the diquark form-factor gives a good description of the FNAL 
data near midrapidity over 2<p t < 10 GeV [7j. Thus, the study of high-p t baryon spectra in 
high-energy p(d)+A collisions is very challenging and provides an opportunity for revealing 
soft successive processes from proton breakup to their parton fragmentations. 

It is also known that physics toward forward rapidity of the incident proton from 
midrapidity is accompanied by a change of particle production mechanism: transition from 
string breaking region, where the particle production is dominated by quark- ant iquark pairs 
created from vacuum, to projectile fragmentation region, where initial charge and isospin 
of the projectile are most likely conserved for the particle production. A signal of such 
a transition appears, for instance, in antiparticle-to-particle ratios from p+p collisions at 
large transverse momentum, which are independent of rapidity at less than y ~ 1.5 but 
above this point decrease with rapidity, irrespective of the detected hadron species jS] . In 
the forward-rapidity region, therefore, diquarks carrying away a large fraction (~ 2/3) of 
the incoming proton momentum will increase the importance of its role for leading baryon 
production in the projectile fragmentation region [5|. 

In very high-energy p+A collisions, in particular for the kinematics realized at forward 
rapidity, a dense target close to the black-body (or unitarity) limit should destroy com- 
pletely the coherence of the projectile partons since all of them experience large transverse 
kicks on the order of the saturation scale (onset of gluon saturation phenomenon) and hence 
fragment independently, mainly into pions jTU] • This leads to a strong suppression of the 
forward baryon number. Also, according to this idea the p t spectrum is expected to be 
rather flat up to the saturation scale. Although recent BRAHMS d+Au data at y ~ 3 [3] 
are almost flat over the range 1 < pt£ 3 GeV, these data do not seem to support this idea 
of forward baryon number suppression at RHIC energy. Rather, the ratio p/ir + ~ 0.8 is 
substantially larger than for p+p collisions x in the range 1 < pt ^ 3 GeV. Therefore, addi- 
tional quantitative theoretical studies of forward baryon production in p(d)+A collisions, 

1 Pythia simulations ^J also underpredict this ratio by a factor 2 ~ 3 as compared to the data [3]. 



and its relation to saturation physics is of importance in detailed comparisons with the 
data. 

Another important experimental observation is that the ratio of hadron ^-distributions 
in d+Au versus p+p collisions is suppressed at large rapidity [12] , while showing a (species- 
dependent) Cronin enhancement at midrapidity [T3*lll4| . Within the Color Glass Conden- 
sate (CGC) formalism [TS], the disappearance of the Cronin peak at moderate p t is induced 
by quantum evolution of the dense coherent gluons in the target [THl H3 [TH] , therefore sig- 
naling the emergence of saturation physics 2 . In the CGC picture, the overall growth of the 
gluon density is decelerated with rapidity because the rate of gluon fusion becomes com- 
parable with that of gluon emission. With increasing rapidity, thus the cross section of the 
gold target at fixed impact parameter grows less rapidly than that of a proton target [21] . 

Along the lines of the CGC formalism, Ref. [22] developed an asymmetric 
DGLAPp ro j [2H] <8> BFKLtarg [21] factorization scheme which accounts for recoil effects due 
to collinear gluon radiation in the projectile, because radiation becomes important with 
higher p t |25j. Furthermore, in Ref. J2H] a parameterization of the anomalous dimension 
7 describing the quantum x-evolution away from the Glauber-Mueller [2Z] or McLerran- 
Venugopalan saturation models [2H] (which do not include evolution due to fixed 7 = 1) 
over a wide range of rapidity was provided. Indeed, these formulations give a good descrip- 
tion of charged hadron or neutral pion production. Then, the importance of recoil effects 
and the appropriate use of the p t and ^/-dependent anomalous dimension in the high p t 
region were discussed quantitatively [221 I2H] 

In this paper, we discuss vr* 1 , p, p production with high transverse momentum (p tr > 1 
GeV) in the deuteron fragmentation region (y = 3.0) of d+Au collisions at RHIC energy. 
At such rapidity the light-cone momentum fractions of the partons participating in the 
scattering becomes very asymmetric between projectile (d) and target (Au) [2H]; large- x 
quarks with x v = O(10 _1 ) from the deuteron collide with (many) small- x gluons, xa = 
O(10~ 3 ~ 10~ 4 ) from the gold target (see figures in |2*2ll2fci| ). Therefore, one should treat 
the deuteron as a dilute and the gold nucleus as a dense object involving saturation effects. 
Its saturation scale is given by Q 2 . ~ A l ' z e Xy with A ~ 200 the mass number and y the 
rapidity of the target gluons. A is a constant describing the growth rate of the saturation 
momentum with 1/x. Thus, the saturation effects in the target are most evident at large 
y and the black-body limit of the target discussed in JTU] is most likely expected there. 
On the other hand, hadron production close to the deuteron beam rapidity is dominated 
by its valence quarks (or their diquark bound state). Energy loss of these leading quarks 
traversing the cold nucleus might occur beyond the leading log Q 2 approximation J3U] ■ 

We will basically use the same expressions for the single-inclusive hadron production 
cross sections as discussed in Refs. [221 EH] in this paper, but will reformulate them slightly 
to include mass corrections (massless particles have been assumed in Refs. [2*21 l2o]); cf. 
Sec. [21 The collinearly factorized form derived here represents our basic equation. In 
subsequent sections we discuss details of the three "building blocks" (i.e. parton distribution 
function (PDF), fragmentation function (FF) and dipole cross section), which are universal 

2 Other scenarios also were proposed to explain the phenomenon [T§10U| . 



(or process-independent) objects. As emphasised in [7j, the scattered constituent diquarks 
represent a possible source of the baryon excess measured at semi-hard p t . This forces us 
to incorporate the diquark component into the first building block, i.e. into the PDF of a 
proton. We shall employ the existing diquark PDF from Ref. [31], which is based on the 
extreme (scalar) diquark picture with Q 2 -dependent form-factor of the diquark, to discuss 
the p t spectra of hadron cross sections including scaling violation. Here we will introduce 
the breakup probability of the diquark scattered off the CGC assuming that the cross 
section for diquark scattering is damped by a Q 2 -dependent form- factor |7j. Regarding 
diquark fragmentation into hadrons, phenomenological considerations using counting rules 
have been applied to the data, and also more theoretical models without scaling violations 
have been proposed [321 [7|. We will, however, not rely on these models but rather model the 
diquark FFs in terms of the standard FFs with scaling violations. As our standard baryon 
FFs we employ parameterizations extracted from JETSET simulations [33], rescaling the 
Q 2 -dependence such as to obtain a good fit to the KKP [31] or AKK 35J FFs. 

Although considerable model-dependences originate from mainly the scaling violation 
of baryon or diquark FFs, the inclusion of constituent diquarks at large x in the deuteron 
proves essential to explain the measured baryon excess relative to standard PDF and FF 
sets. For tt^ and p production, on the other hand, the diquark contribution is not so large. 

This paper is organized as follows. In the next section, we formulate the single-inclusive 
hadron production cross sections in high-energy d+Au collisions including mass corrections. 
In Sees. [3J and HI the predictions from the the model are confronted with experimental data 
for forward ^ and p/p production, respectively. A detailed discussion is given in Sec. El 
Finally, we summarize in Sec. [HI 

2 Forward hadron production in proton-nucleus col- 
lision within the CGC formalism 

It is well known that collinear factorization theorem is generally broken by higher-twist 
(multiple scattering) effects in hadron-hadron collisions, although the validity of the the- 
orem is proved up to the twist-4 level J3H]. Within the CGC formalism, however, all 
higher-twist effects can be included into one building block, i.e. the dipole forward scat- 
tering amplitude. Hence, even if one considers higher-twist effects in p+A collisions, the 
predictive power of the theorem is still preserved [37]. It was verified also for the case 
with one-loop radiative corrections to projectile partons [22]- This expression [22] is valid 
only under the condition that the target is in the saturation or so-called "extended geo- 
metric scaling" regime while the projectile is dilute [38]. Therefore, this approach does 
hold for forward rapidity kinematics which we are interested in here, unlike the usual k t - 
factorization approach [35]. In this section, we shall briefly derive collinearly factorized 
forms describing the single-inclusive hadron production cross section in p+A collisions. 

To describe the p + A — > h + X process depicted in Fig. Q we work with light-cone 
momenta in a frame where the proton has large plus component P p = (ws/2, 0,0*) and 



p 




Figure 1: Kinematics of the p+A —>■ h+X process. The incoming parton i from the proton 
interacts with the target CGC field of small- £ gluons to all orders via a 2 — > 1 process. 
The outgoing parton / hadronizes into a massive hadron h with momentum fraction Zh- 

the nucleus a large minus component Pa = (0, Js/2,0 t ). y^s denotes the center of mass 
energy. The incoming parton % carries a fraction x p of the proton momentum P+ in- 
teracting with the incoming small-x gluons of the dense nucleus (the CGC) through a 
2 — ► 1 process. After the eikonal scattering process, the outgoing parton / has mo- 
mentum pf = (q t e yf /\^2,q t e~ yf /^/2,q t ) with rapidity yj and produces a massive hadron 
h with the momentum fraction Zh and mass m^. The detected hadron has a momen- 



tum ph = (m t e yh /\^2,m t e~ yh /^/2,pt), where m t is the transverse mass m t = Jml + Pt 
and Dh is the hadron rapidity. The Feynman-x of the produced hadron is given by 
xp = Phi Pp — m t e Vh J ' y/s. In the eikonal approximation we have Zh = Vhlv) = X p/ X F be- 
cause of x p = q t e ys I ' s/s. The internal variables of the parton /, i.e. qt and yf, therefore, can 
be expressed in terms of p t , m t , y h and x p : q t = p t /z h = p t x p /x F and y f = y h + log(m t /p t ). 
Another important internal variable is the momentum fraction xa carried by the small-x 
gluons in the nucleus. It is related to that of the impinging projectile parton x p in the 
eikonal approximation by xa = x p e~ 2yf (see Appendix B in [22] )• This leads to the rapidity 
of the gluons, 

y A = log(l/xA) = Iog(l/x p ) + 2y h + \og(m 2 t /p 2 t ). (1) 

Here, complying with standard practice, we have defined the direction of the final hadron 
produced in the "proton fragmentation region" as the positive ^-direction, i.e. yh > 0. 

In the framework of collinear factorization, the single-inclusive hadron production cross 
section can be written as a convolution of the inclusive parton / production cross section 
with the PDF of the initial parton i, fi/ p , and with the FF of / into a hadron h, Dh/f'- 



da(pA -> hX) _ da(pA -> hX) 
7rdyhdp 2 d 2 b -ndy h dm 2 d 2 b 



da(iA -> fX) 



= J2J dx p dz h dy f dq 2 f i/p (x p , Q 2 ) — - 2 ^ D h/f (z h , /ij) 

x 5{m 2 - M(z h , q t ) 2 ) 5{y h - Y(y f , z h , q t )), (2) 



where M(z h , q t ) 2 = m 2 h + (z h q t ) 2 and Y(y f , z h , q t ) = y f - \og(m t /z h q t ) as derived above, b 
is the impact parameter with respect to the center of the nucleus. The factorization scale 
of the PDF is denoted by Qf and the fragmentation scale of the FF by [if, both of which 
are hereafter set to Qf = \if = p t . 

After integration over yf, q 2 and Zh, we obtain the impact-parameter averaged single- 
inclusive hadron production cross section, 

dajpA -> hX) = ^j l dx ^ [ R * dbb 
ndyhdp 2 2nS A Jx f p xf Jo 



fq/p(x p , Q f ) N F [ — p t , y A , b 




x 



+ fg/p(x P , Qf) N A [—Pt, Va, b ) D h/g [ — , \if ] (3) 

Here we used da{iA -^ fX)/irdyfdq 2 d 2 b = qfS(qJ — x p Pp)N(q t ,yA, b)/(27i) 2 , where 
N(q t ,yA,b) is the scattering probability of dipoles from the nucleus. Np corresponds to 
a projectile quark impinging with the target small-x gluons while iV^ is for a projectile 
gluon. The indices q in f g / p and D h / q are summed over all quark species. For simplicity 
we assume a spherical nucleus of radius R A with a sharp edge, where its profile function is 



given by T{b) = 2JR\ — b 2 . Sa is the transverse area of the nucleus (= vri?^). 

The same expression as Eq. Q was derived in Ref. [22] assuming massless hadrons 
and hence this equation is consistent with that of [22] i n the massless limit. To leading 
\ogp 2 accuracy it was verified in [22] that the recoil of the projectile parton by hard gluon 
radiation can be converted into the Q 2 -evolution of the PDF and FF according to the full 
DGLAP [23] evolution equations. The recoil effect are important for the interpretation of 
the forward-rapidity data from RHIC. 

The dipole profiles of transverse size r t at an impact parameter b are defined in the 
fundamental and adjoint representations of SU(N C ), respectively, as 

N F (r t ,y A ,b) = ^Tr c (l-V\b-r t /2)V(b + r t /2)), 

N A (f t ,y A ,b) = jJ-^Tr c (l-rt0-<r t /2)U(b + 7 f t /2)), (4) 

where V and U denote Wilson lines along the light cone [40J in the corresponding repre- 
sentation, and their correlators are averaged over the color source in the nucleus. 

For the actual description of the dipole profiles we rely on the KKT model [17J, which 
is a phenomenologically reasonable parameterization facilitating comparisons with exper- 
imental data, rather than on solutions of the intricate JIMWLK equations jUJ], which is 
not yet feasible. For a dipole in the adjoint representation, 



N A (n,yA,b) = 1 -exp 



—Ar 2 t Q 2 s {y A ,b)) 



l{rt,VA,b) 



(5) 



where Q s {y A , b) is the saturation scale of the target nucleus and j(r t , y A , b) the anomalous 
dimension of its gluon distribution with saturation boundary condition. 7 = 1 reduces (0) 



to the Golec-Biernat-Wusthoff saturation model jH] or the classical MV model [28J. The 
dipole in the fundamental representation, N F , differs by a factor of Q 2 — > Q 2 Cf/Ca = 
±n 2 IT7I 

The saturation scale Q s at b = is given for a nucleus of mass number A (Au(197)) 
as [22] 

Q 2 (y A ,b = 0) = A^Ql (£°V = A^Qlxyy-, (6) 

\x A J 

where Q ~ 1 GeV, A ~ 0.3 and x ~ 3.0 x 10" 4 are fixed by the DIS data |H]. 
The energy dependence of Q s is controlled through the constant growth rate A = 
d\og(Q 2 /AQ CD )/d\og(l/x A ), which is obtained from fixed-coupling LO BFKL evolu- 
tion [12] • Since the squared saturation scale has dependence on the impact parameter, 



for instance, through the nuclear profile T(b) = T(b = OWl — (b/R A ) 2 in the hard sphere 
approximation of the nuclear target and the pointlike proton projectile, we define the 
saturation scale as 



QUva, b) = Qi(y A , b = Otyl - (b/R A y. (7) 

This naive approximation for the nuclear surface is sufficient for minimum bias observables. 
The anomalous dimension 7 is parameterized as J2EJ 

/ n uh \ \\og(l/r 2 Q 2 (y A ,b))\ 



\y A +\\og(l/r 2 Q 2 s (y A ,b))\+d Jyl 



where 7 S ~ 0.627 is the anomalous dimension for BFKL evolution [21] with saturation 
boundary condition, i.e. for evolution along the saturation line |42j . and d is a free parame- 
ter which is fitted to experimental data. Throughout this paper, we will make replacement 
j(rt,y A ,b) — > j(l/q t ,yA,b). As indicated in [2E], this parameterization of 7 stays within 
7 = 0.6 ~ 0.8 at large rapidity yh over a comparatively wide range of p t = 1 ~ 5 GeV 3 . 

Below we focus on the minimum-bias cross section obtained by impact-parameter aver- 
aging of Eq. ([2]). Since in the integrand of (J2J) the impact-parameter dependence is carried 
only by the saturation scale, we take the average value of Q 2 (b) with respect to b instead 
of integrating Q over b: 



(Q 2 s (b)) = £. f A db 2 Q 2 s (y A , b) = \ Q 2 (y A , 6 = 0). (9) 

DA Jo 3 



We have used Eq. (J7|) where this approximation is valid to good accuracy 



3 As emphasized in [21], this remains inside the (at least, extended) geometric scaling regime. When 
going beyond the scaling regime like higher-p t or lower rapidities, Eq. Q assumes that crossover from the 
scaling regime to the perturbative one is very slow and smooth. A similar behavior of 7 was discussed in 
Ref. g3j 

4 Refs. [2U EH employ an effective saturation scale for minimum bias collisions such as {Q 2 s {b)) = 
A e ffQo( x o/xA) X with A e ff = 18.5, where the factor AA* is by ~ 30% smaller than the corresponding 
factor (2/3)A 1 / 3 in ©. 



Then, Eq. (J3J) reads 
da m . h {pA -> hX) 1 



ndy h dp 



(2k) 



■'-F 



dx p ^- 

X F 



f q/p {x pi Q))N F (^p t ,y A , (Q 2 s (b))) D h/q [?*.,$ 

\ X p / V Xr) j 



+ fg/pfo, Q))N A (%,, y A , (Ql(b)}) D h/g fe n) 

\X p / \ Xr) , 



v, Jb P / \ JU r) 

where the Fourier transform of the dipole profile functions is given by 
N A , F (qt,y A ,(Q 2 M) = - fd 2 r t e l ^N AF (r t} y A ,(Q 2 s (b)}) 



-2tt / dr t r t J (r t q t ) N A>F (r t , y A , (Q s (b))). 
Jo 



(10) 



:n^ 



3 Forward 7r ± production in d+Au collision at RHIC 

We are now in the position of applying our results (jlOJI to deuteron-gold collisions at RHIC 
energy yfs = 200 GeV, at large rapidity (y^ = 3.0). We first address minimum bias 
tt^ production observed by BRAHMS collaboration 03]. This process is dominated by 
valence quarks since x p > 0.1 on the deuteron side and hence might be affected by initial 
quark correlation (like diquarks) with large x momentum inside deuteron. 

In this section, our interest is twofold: First, we check the validity of our leading order 
CGC formalism using standard parameterizations of the PDF and the FF, and determine 
the free parameter d as well as the if -factor prescribed to NLO corrections to reproduce 
well the experimental data. Second, using those very same parameters (d and K) we 
investigate the diquark contribution. 

For a deuteron projectile, we treat its parton distributions as a simple superposition 
of those in a proton and a neutron without any nuclear modification, which indeed has 
negligible effects (less than 5%) for the deuteron. The parton distributions in the neutron 
are obtained through isospin symmetry, f U) u/ p = fd y d/n an d fd,d/ P = fu,u/n, and other parton 
species have the same distribution as those in the proton jlHl • This isospin symmetry may 
affect the final pion productions composed of flavor SU(2) light quarks, where it leads to 
no difference between light 7i + (ud) and n~(ud) productions to good accuracy. Actually, 
such behavior can be seen for p t < 4 GeV in the BRAHMS data plotted in Fig. El 

We first use the leading order (LO) CTEQ5 PDFs jUj. Then, we have several choices 
for the FFs, e.g. the Kretzer gg], KKP 34 and AKK ^ sets. The FFs by Kretzer 
assume that the charge-conjugation symmetry D q / W ± = D q /^ holds at low input scale. 
Both KKP and AKK FFs (the latter provide flavor- dependent FFs of light quarks as an 
update of KKP) provide only the average of charged pion FFs. Therefore, all three FFs 
give the same neutral pion FF and should show almost the same p t - spectra at LO. 

We checked that the result with the LO Kretzer FFs gives a very similar curve to that 
of the LO KKP FFs (but different if-factor) over the range 1 < p t < 5 GeV, with a pion 
mass of m n = 0.14 GeV. The result with the LO Kretzer FFs is in good agreement with 
the data at d = 0.6 and K = 1.4, while the LO KKP FFs require d = 0.6 and K = 1.0. 
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Figure 2: p t spectra of tt + and ir~ from d+Au collisions compared to BRAHMS preliminary 
minimum bias data (yh = 3.0). The lines show the CGC results with the LO Kretzer FFs 
and K = 1.4, where either the LO CTEQ5 or diquark PDFs are used. 



The main difference is due to the way that gluons fragment to pions; its contribution in 
the KKP set is larger than that in the Kretzer one [2] . The sensitivity of the results to the 
value d is not so large in the range of d = 0.6 ~ 1.2 and is visible only at high p t , where 
the data however have large errors. In Fig. |2]we show the LO Kretzer result with K — 1.4 
as the dotted line. 

Next, in order to investigate the diquark contribution to this process, we make use of 
the PDFs given by the Stockholm diquark model [SI] . It assumes that only scalar diquarks 
are genuine bound states and other axial-vector diquarks are negligible. As is well known, 
the diquark is classified into 3 C and 6 C channels in color SU(3). The scalar diquark belongs 
to the 3 C channel while the axial- vector diquark is in the 6 C representation. The reason why 
we adopt only the former is that the one-gluon exchange between quarks is attractive in the 
3 C channel while repulsive in the 6 C channel. Thus, for a first estimate it is rather plausible 
to focus on the more tightly bound scalar-diquark. This parameterization is obtained from 
a fit to the proton structure function F% observed in high energy e+p collisions at SLAC, 
BCDMS and EMC, over the wide range 1 < Q 2 < 200 GeV 2 . 

Ref. [3T] has two possible parameterizations which differ by the choice of a mass scale 
in the diquark form-factor. We choose the parameter set with M 2 = 10 GeV 2 and a dipole 



form-factor 

F2{Q '> = (n-qSw (12) 

It ensures that with larger Q\ all diquarks are resolved and the usual quark PDFs are 
recovered. The scale M is related to the internal binding energy of the diquark and 
M 2 = 10 GeV 2 corresponds to a diquark radius of about 0.2 fm, which is obtained from the 
mean-square radius (r 2 ) = 6dF 2 (Q 2 )/dQ 2 \q2 =0 . It yields larger contribution from diquark 

component due to its tightly bound state than that of M 2 = 3 GeV 2 . The parameterization 
for each parton, then, is given as [3*T] 



XpUyXp, vjj 



x 



(0.52 + 0.37s)x p A5+0A2s {l - x p ) 2 - 2+2 - 39s - x° p 71 -°- 02s {l - x p ) 7 - 3 -°- 51s 
(0.52 + 0.375)5(0.45 + 0.12s, 3.2 + 2.39s) - £(0.71 - 0.02s, 8.3 - 0.51s) 

™0.9~0.83sfi _ „ \5.0-1.88s 

, JZ2 \z p) \1-F 2 (0 2 )] 

£(0.9 -0.83s, 6.0 -1.88s) l Wf,b 

0.9-0.83s/-| _ N5.0-l.88s 

XpSyXp, Uc j ) = XpUyXp, V/J = XpClyXp, VjJ = %pS\%pi Usf) 
= (0.35- 0.06s) (1 - x p ) 6 - 89+a75s , 

2,0.93-0.523 /•-]_ _ x \1.5-l.ls 

/ ' ,i " M ' / ' ^ } :::: £(0.93- 0.52s, 2 ? 5- 1.1s) * ' Wfh 



Jd Q {x Pi Q 2 s ) = D ; nnQ "" TT^ F 2 (Q% (13) 



where we assumed flavor 577(3) symmetry for the sea quark distributions, and 
Jdq denotes the diquark distribution, B(ji,v) Euler's Beta function and s = 
log[log(Q^/A 2 )/log(Q 2 /A 2 )] with A = 0.2 GeV and Q 2 = 4 GeV 2 . 

When one views the d+Au collisions in this diquark picture, one faces other theoretical 
uncertainties besides modelling of the diquark distribution in the deuteron: the breakup of 
diquarks scattered off the CGC and the diquark fragmentation into hadrons. The former 
depends on the collision dynamics between diquarks with a finite size and the CGC. To 
take into account this breakup probability in a simple way, we express it as the dipole 
form-factor of (|12j) with M 2 = 10 GeV 2 as done in Ref. [7j. This probability goes to zero in 
the high-p 4 limit (Qf = p t ), and then the diquark bound state is completely broken inside 
the CGC, where the constituent u and d quarks fragment independently. 

Although so far various models for the diquark fragmentations have been proposed [32] , 
they have focused mainly on how the FFs behave as a function of x, leaving aside scaling 
violation by the Q 2 evolution. For our present investigation of p r spectra of pions and 
baryons from diquarks it is important to account for the scaling violation. For this purpose, 
as a naive model estimate we express the diquark FFs in terms of the LO Kretzer FFs, 
which describe fragmentation into charged hadrons separately 5 . We approximate the 

5 The standard FFs like the Kretzer set describes fragmentations from a color triplet parton. Since the 
scalar diquark is in the color anti-triplet representation, it may hadronize similarly to a parton, at least as 
far as color is concerned, if the size of diquark is ignored. 
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(-U(i)-diquark FFs into charged pions as 

Dn+(ud)/DQ(ud)\ z h, V>f) = D^/d { z h,^f)F (fJ,f) 

+ [D«l; u L °{z h , /*?) + D^/ d LO (z h , ^)](1 - F 2 (/iJ)), 

Dn-(ud)/DQ(ud)(Zh, Hf) = D w lj u (z h ,fI f )F (/J, f ) 

+ [D£jf>(z h ,fi) + D«IJ u LO {z h y f )){l - F 2 (/iJ)), (14) 

where we multiplied the FFs by the breakup probability. The first term is the contribution 
where the valence quarks of the produced pions do not include neither one of the quarks 
from the scattered diquark, and thus the pion production occurs via quark-antiquark pair 
creation from the vacuum, like D n +/ d or D n -/ U . The second contribution corresponds to 
diquark breakup followed by independent fragmentation of each quark into pions. 

In Fig. |21 we plot the diquark result, using Kretzer FFs both for the diquarks as well 
as for the single partons, as the dashed line (with d = 0.6 and K = 1.4). This curve, 
however, has a problem, since it should approach the result with CTEQ5+Kretzer (dashed 
line) at high p t due to F 2 (p 2 ) — > 0; that is, the diquark contribution should vanish in this 
limit and the single parton picture should be restored again, but apparently the former 
is much harder than the latter. To improve this unfavorable behavior, we simply change 
the dipole form-factor F 2 (p 2 ) to F A (pf) for both the PDFs and the breakup probabilities. 
This result is shown as the solid line, which is softer at high p t and is in good agreement 
with the CTEQ5+Kretzer result, especially at high p t . Since this change of the form-factor 
increases the diquark radius from 0.2 fm to 0.3 fm 6 , the new form-factor becomes more 
sensitive to pt and the diquarks are broken up more easily. The solid line is still somewhat 
harder than the dotted one at low p t , but it remains a good fit to the data. Thus, the 
result in the diquark picture can reproduce well the data within error bars. As expected, 
this implies that the diquark contribution to pion production is not very large. 

4 Forward p and p production 

In this section we examine forward proton or anti-proton production in d+Au collisions at 
RHIC and apply Eq. (jlOj) to obtain the cross sections within the same setup as for pion 
production (discussed above). At the RHIC energy y/s = 200 GeV and in the deuteron 
fragmentation region (y^ = 3.0), we first show the BRAHMS preliminary data [12] for 
minimum bias cross sections of p (solid circle) and p (open circle) production in Fig. El 

Unlike for pions, charge symmetry is broken more strongly for p, p production. The 
proton cross section exceeds that of anti-protons by one order of magnitude for l.5<pt£3.5 
GeV. This feature can also be seen in p+p collisions at large rapidity and high p t [15]. Such 
a large asymmetry between them, and in particular the high abundance of protons, may 
be associated with the production mechanism involving the fragmentation of diquarks into 
protons. 

6 This radius is still smaller than the other parameter set of [7] with M 2 = 3 GeV 2 , where \J {r 2 ) ~ 0.4 
fm. 
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Figure 3: p t spectra of p and p in d+Au collisions from BRAHMS preliminary minimum 
bias data (yh = 3.0) compared to computations within the CGC formalism (K = 1.0). The 
lines show the results with the common LO CTEQ5 PDFs and the Lund FFs rescaled to 
either the KKP or AKK FFs with K = 1.0. The lines with symbols show the results with 
the diquark PDFs. 

To our knowledge, the best parameterization of baryon FFs handling the Q 2 evolution 
to some extent is given in Ref. |33j . where the FFs are obtained by parameterizing results 
computed by Monte Carlo simulations in the Lund string model, evaluated at invariant 
mass W = 2Qf. The advantage of using this parameterizations is that we can treat pro- 
ton and anti-proton production separately while the KKP and AKK FFs only provide 
averages. Moreover, the Lund string model JETSET also models mass effects in the frag- 
mentation which are important when the momentum of the outgoing parton is not much 
larger than the mass of its daughter hadrons. However, the Q 2 r dependence of these Lund 
FFs D B / q (zh,Q'f) has not been implemented very carefully: Ref. J33J only requested that 
the average multiplicity Jq dzhD B / q (z hl Q 2 ) satisfies the Q 2 dependence of the experimen- 
tal data. This treatment would not describe correctly scaling violations, especially in the 
high-Zh regions, where the FFs change steeply. 

To check the Q 2 t behavior, we computed the average pt-distribution of protons and anti- 
protons, (p + p)/2, and compared with the results obtained with the NLO KKP and NLO 
AKK FFs, using the LO CTEQ5 PDF for all three FF sets 7 . The ratios KKP-NLO/Lund 



7 The reason why we employ the NLO KKP or NLO AKK FFs is as follows: Recent STAR data show 



12 



and AKK-NLO/Lund are shown in Fig. 0J where the baryon mass is set to m p p = 0.935 
GeV and we choose d = 0.6 and K = 1.0. 



re 3 
cc 



1 - 



(CTEQ5 + KKP-NLO)/(CTEQ5 + Lund) 
(CTEQ5 + AKK-NLO)/(CTEQ5 + Lund) 




1 1.5 



Figure 4: Summed p + p yields obtained with the NLO KKP or NLO AKK FFs divided 
by the Lund FFs. 

We see that the ratios change steeply as a function of Qf = p t and that the Q"j depen- 
dence of the Lund FFs is significantly harder than those of KKP and AKK sets. In more 
detail, the cross section in the Lund model is larger than that with the KKP set for p t >, 1.5 
GeV, while the behavior is opposite for the AKK set at p t ^ 3.5 GeV. This forces us to 
improve the Lund model FFs to enforce a Q ^-dependence that matches that of the KKP 
or AKK FFs. In this figure we can also see that the KKP FFs are smaller by a factor of 
5 ~ 12 than the AKK FFs in the region 1.5<,p t ~ 5 GeV and the deviation becomes larger 
at higher p t . 

For that purpose, we fit the ratios using each function via 



KKP(NLO)/Lund 
AKK(NLO)/Lund 



-O.Ol + Got/1.56) 1 - 36 ] 2 - 01 ' 
1 



[0.12 + (pt/3.70) 1 



9111.38 ' 



(15) 



that high-pt p + p yields around midrapidity in d+Au collisions at y/s = 200 GeV agree better with NLO 
pQCD calculations using the AKK rather than the KKP FFs |2|. The latter is known to underestimate 
the data by one order of magnitude and so it would be best to adopt the AKK FFs for baryon production. 
Unfortunately, the AKK FFs only provide a NLO parameterization. The use is inconsistent within our LO 
formalism, but despite such inconsistency, the use of AKK FFs for baryon productions is of importance 
to reproduce the yields. We also checked that the NLO KKP FFs are quite similar to the LO KKP FFs, 
but the former is a bit softer than the latter at high p t . 
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These functions provide rather reasonable fits. Multiplying the Lund FFs by the fitting 
functions (J15JI . we then rescale the Q 2 r dependence of the FFs. We shall employ the same 
rescaling functions for both p and p FFs. 

For the diquark fragmentation into baryons we construct the following model in the 
same way as for pion production: 

D p ( uud )/DQ(ud)(Zh, Vf) = D-n+iudb/d^hi ^f) F W) 

+ [iffiizH, fi) + D^(z h , /4)](1 - F 2 ^)), 



Dp(uud)/DQ(ud)( z h,Hf) — TvPffu + Dp/T W (Vf 



Lund I _, ,,2\ i r\Lund( _, ,,2\l/i ZT^/,^' 



+ [Iffi(z h ,ti) + Iffi{z h ,i4)](l-F'Q4)), (16) 

where D Kr - LO and D Lund denote the LO Kretzer and Lund sets respectively. The first 
terms express the diquark fragmentation as a single entity and the second ones do as a 
system consisting of independent valence quarks. The first term of the proton FF is based 
on the conjecture that the fragmentation of the (w<i)-diquark proceeds via pick-up of a u- 
quark from the vacuum, similar to that of d into 7r + , if the finite size of diquark is neglected. 
The first term of the anti-proton production treats a diquark as a single entity (like u or d 
quarks) and hence takes their average. 

In Fig. El we plot the results of Eq. (|1UJ1 rescaled to the KKP or AKK FFs, and using 
the form- factor F A . For anti-protons, we show a total of four curves with the CTEQ5 or 
diquark PDFs, rescaled to either the KKP (dashed) or AKK (dotted) FFs. As indicated 
above, the KKP set significantly underestimates the anti-proton data by about a factor 10, 
independent of the choice of PDFs. The AKK FFs fit better although the data appears 
somewhat softer than our results. The difference between the CTEQ5 and diquark PDFs 
is generically small but increases at low p t , similarly to the pion case discussed in Sec. El 
Both these curves are in good agreement with the data within the large error-bars. 

For protons, we first show the CTEQ5 result rescaled to the AKK set as the solid line, 
which underestimates the data by a factor 2 ~ 3 in the low-pt region, but approaches the 
data at high p t . Next, using the same FFs, the diquark result is plotted as the square+solid 
line, which enhances the cross section at low p t , but still underestimates the data. This 
feature is, however, unlike the cases of anti-protons or pions and the enhancement of the 
cross section at low pt is favored by the data. As a somewhat extreme case we also plot a 
curve corresponding to "complete diquark survival" by setting the breakup probability to 
zero (i.e. F A (fi 2 r) = 1) as the open-circle+solid line. This result is very close to the data in 
the low-pt region. 

The p/tc + and p/tc~ ratios will be relevant for investigating the pt dependence of baryon 
production more precisely. In Fig. |3] we plot the data extracted by taking the ratios of 
the corresponding data points from Figs. El Here, as seen in Fig. |21 there are no data 
points of the ir~ in the range 1.5<p t ^, 2.1 GeV. To interpolate the lack in this interval (6 
points), we substitute the corresponding data points of the 7r + . This is reasonable because 
for a deuteron projectile isospin symmetry seems to work well at low pt as discussed above. 
This figure shows that the p/ir + ratio is almost 0.6 ~ 0.9 which is almost flat, while the 
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Figure 5: Ratios of p/^ + and p/vr~. The experimental points were obtained by taking 
the ratios of the corresponding data from Figs. El and El The lines show the results with 
the common LO CTEQ5 PDFs for both baryons and pions, where the former uses the 
Lund FFs rescaled to the AKK FFs with K — 1.0 and the latter the LO Kretzer FFs with 
K = 1.4. The lines with symbols show the results with the diquark PDFs. 

p/n~ ratio is 0.04 ~ 0.07, smaller by one order of magnitude, where at p t < 2.5 GeV we 
see no significant ^-dependence, but above this point the ^-distribution starts to decrease 
largely with p t . 

For the p/ir~ ratio, we plot two results: The dotted line uses the common LO CTEQ5 
PDFs for both p and 7r~, and the Lund FFs rescaled to the AKK set for p and the Kretzer 
FFs for 7r~, with different i^-factors: K = 1.0 for p and K = 1.4 for 7T~. The cross- 
dotted line replaces the PDFs by the diquark set, with otherwise the same "ingredients". 
These curves explain the data quantitatively, especially at low p t , with small discrepancies 
between the CTEQ5 and diquark PDFs. The monotonic decrease with p t is supported by 
the data except the last data point. 

For the p/ir + ratio, we show the same curves 8 . First, we employ CTEQ5 PDFs with 
the Lund FFs rescaled to the AKK set for p and to the Kretzer FFs for 7r + (solid line). 
Second, we switch to the diquark PDFs with the same sets of FFs. The experimental data 
cannot be explained by the use of the CTEQ5 PDFs any more, which is too low by a factor 



8 This ratio is above unity around pt = 1.5 GeV, because of the relatively large underestimate of our 
7r + cross section as compared to the tt + data. 
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of 2 ~ 3 at low p t , as already seen in Fig. |H1 The result is almost flat over a wide range of 
Pt and seems to work well only in the high-p t region. Switching to the diquark PDFs may 
explain at least partly the large baryon excess seen in the data at low p t . Diquark effects 
disappear quite rapidly with increasing p t which is not really supported by the data. The 
prediction of the "surviving diquark" picture is moderately larger at low p t and perhaps 
offers a reasonable explanation of the data. 

5 Discussion 

The physical process of d+Au collisions illustrated in Fig. ^ shows that the cross section 
is factorized into three parts, i.e. into a convolution of the form f q>g /d <8> N F) a ® ^h/q,g like 
in Eq. (jlUj) (refer to ref. [22] for verifying this factorization to leading logarithmic accu- 
racy 9 ). These three functions are universal, independent of the processes and associated 
with infrared hadron or nuclear structures, for example with respect to the fact that we 
employ the same forms of Np t A for meson (D M / qg ) and baryon (D B / q>g ) production. In 
other words, the p t spectra may be determined by universal dipole profiles, irrespective 
of detected hadron species. These dipole profiles Np,A describe the interaction of hard 
partons scattering off the fields of nucleus. Bremsstrahlung occurring either "before" or 
"after" the propagation through the nucleus yields logarithmic radiative corrections, which 
are absorbed into the DGLAP (^-evolution of f q , g /d and D h / W |22j. 

As seen from Figs. 121 El the data displays a common power-law behavior over the range 
of l£jpt£j5 GeV, not an exponential form predicted by the parton recombination model [TU] . 
In more detail, it is close to the exponential until pt ~ 2 — 3 GeV, but above this it appears 
rather close to a power-law. In the CGC formalism, the dipole profile displays a power-law 
behavior for the ^-distribution at high rapidity, whose specific shape is governed by the 
anomalous dimension 7 (221 EHj • 

In fact, our dipole profiles give even slightly harder ^-distributions than the data for 
both p and p if we employ the CTEQ PDFs. If the factorization and the CGC framework 
work well, however, such a deviation from the data will be due mainly to the Q 2 -evolution 
FFs (here the Lund baryon FFs rescaled by the AKK FFs). If this is the case, we should 
determine more carefully the scaling violation of the baryon FFs in the Lund string scheme. 
Work on new parameterizations by means of Lund JETSET simulation will be reported 
elsewhere |oT] . 

For p production, a remarkable point is that the diquark picture does have a noticeable 
effect on the spectra at large rapidity. In fact, for the specific parameterizaiton employed 
here, the slope is even steeper than seen in the data. This comes from mainly the p t - 
dependence in the form-factors, which are used both in the PDFs and in the breakup 
probability of diquarks. As a reference, we considered the case of completely pointlike 
diquarks, which are not broken by scattering off the CGC at all. At a first glance this result 

9 The process violating the factorization like the second rescattering of the hard parton with the target 
after the first rescattering and the succeeding gluon radiation vanishes in the high energy limit |49II5(J| . at 
least if we calculate in light-cone gauge. 
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is closer to the data compared to the results with breakup. However, this should be taken 
with caution, because this "surviving diquark" picture maximizes the contribution from 
diquarks and is inconsistent with the baryon suppression scenario in the deep saturation 
regime [TU] . 

For a more quantitative comparison to experimental data we require more precise infor- 
mation of the Q 2 -dependence of diquark distribution and fragmentation functions (which 
can not be measured in e + e~ annihilation like ordinary quark FFs). It might also be possi- 
ble that the deuteron is composed of three-diquark system |52j . In the present formulation 
the third constituent diquark has not been considered. This contribution may lead to 
additional baryon formation. 

6 Summary 

We have focused on forward p and p production in d+Au collisions using universal dipole 
profiles developed in [2211211 within the CGC formalism. This formalism is in a good agree- 
ment with the high p t -spectrum of forward pions at y = 3.0, where its power-law behavior 
is described well by the dipole profile in the formalism. However, the abundance of protons 
can not be explained by the usual PDFs (like CTEQ5) and Lund FFs (rescaled by AKK 
FFs) , while that of anti-protons is consistent with our results. Similarly, large deviations of 
the p/ir + ratio from Pythia simulations is already observed for forward p+p collisions [13] . 
One possible explanation for this extra production of protons at p t — 1 ~ 3 GeV compared 
to standard PDFs is the direct formation of protons from (dominantly scalar) diquarks, 
whose contribution is incorporated into our formulation, including a diquark form-factor 
and fragmentation into hadrons. This diquark contribution is essential only for proton 
production but is less important for pions and anti-protons. 
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